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The creation and manipulation of multipartite entangled states is important for advance- 
ments in quantum computation 1 and communication 2 4 , and for testing our fundamental 
understanding of quantum mechanics 5 and precision measurements 6 . Multipartite entangle- 
ment has been achieved by use of various forms of quantum bits (qubits), such as trapped 
ion£l5|, photons 9 , and atoms passing through microwave cavitie^. Quantum systems based 
on superconducting circuits have been used to control pair- wise interactions of qubits, either 
directljPES through a quantum bu&2E3, or via controllable coupling 16 . Here, we describe 
the first demonstration of coherent interactions of three directly coupled superconducting 
quantum systems, two phase qubits and a resonant cavity. We introduce a simple Bloch- 
sphere-like representation to help one visualize the unitary evolution of this tripartite system 
as it shares a single microwave photon. With careful control and timing of the initial condi- 
tions, this leads to a protocol for creating a rich variety of entangled states. Experimentally, 
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we provide evidence for the deterministic evolution from a simple product state, through a 
tripartite W-state, into a bipartite Bell-state. These experiments are another step towards de- 
terministically generating multipartite entanglement in superconducting systems with more 
than two qubits. 

With the development of quantum information science 1 , entanglement of multi-particle systems 
has become a resource for a new information technology. In particular, three-particle or tripartite 
entanglement allows for teleportation 2 , secret sharing 4 , and dense coding 17 , with connections to 
cosmologjff^. Over the last decade, the development of exquisite control over quantum systems 
has led to various demonstrations of tripartite entanglement^. Genuine tripartite entanglement 
is delineated by two inequivalent classes of stated GHZ (Greenberger-Horne-Zeilinger) and W, 
where the W-state involves only a single photon shared amongst three systems. Utilizing multipar- 
tite entanglement in a solid- state-qubit system has only recently received theoretical attention^^. 
Thus far in superconducting systems, bipartite entanglement has been verified by two-qubit quan- 
tum state tomograph^ 13 and used to perform a quantum algorithm 15 . Spectroscopic evidence for 
three-particle entanglement was observed for two current-biased phase qubits coupled to a lumped 
element LC-cavity as well as for Transmon qubits JHH11 In the experiments described below, we 
first verify the spectroscopic signature of three coupled systems. Next, we demonstrate coherent 
interactions. Frequency detuning of the third system is used to verify the proper change in the 
time evolution of two versus three coupled systems. Finally, we describe a free-evolution process 
and a visualization technique as a method for deterministically preparing arbitrary single-photon 
tripartite entangled states. We present evidence for the proper operation of this protocol by mea- 
suring the time-dependent behavior of the two phase qubits. Here, entanglement is not verified 
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directly, but the data are consistent with theoretical predictions. Proper execution of this protocol 
can prepare the system in a Bell or W-state, as well as arbitrary entangled states. 

In Fig.[T^, we show an optical micrograph of two qubits, qubit 1 and qubit 2, capacitively coupled 
to either end of an open-ended coplanar waveguide cavity whose half-wave resonant mode fre- 
quency is uj c /2h « 8.9 GHz. These cavities have shown coherent properties at the single-photon 
levelpl Flux-biased phase qubitP^can be thought of as anharmonic LC-oscillators in which a sin- 
gle Josephson junction provides enough nonlinearity to address the two lowest oscillatory phase 
states \g) and |e). The energy level separation hujj = E e — E g can be independently tuned over 
a range ~ 7 GHz— 10 GHz on the j-th qubit by use of inductively coupled flux bias coils. An 
additional coil allows us to apply microwave pulses and fast bias shifts, also used for single-shot 
state measurement.^ Independent state readout on the j-th qubit is accomplished by use of an 
inductively coupled dc superconducting quantum interference device (SQUID). We describe this 
system using a two-qubit Jaynes-Cummings or Tavis-Cummings model. 26 In a frame rotating at an 
reference frequency uj^, we approximate the Hamiltonian of the system as 

H = HA c (Ja + [ hA 3 a t a 7 + ih 9j(rf a ~ a](J j)] > 

where the mode operators af and refer to the qubits and the cavity, respectively, with corre- 
sponding detunings Aj/h = ujj — uj^ and A c /h = uo c — uj^. Capacitive coupling C c between the 
qubits and the cavity results in an effective coupling frequency of 2gj/2n « uj c /2ti C c /\/CCj « 
2g/2n ~ 90 MHz for both qubits. The system exhibits decay rates of 71 /2tt ~ 7 MHz, 72 /2tt ~ 10 MHz, 
and k/2tt ~ 1 MHz for each qubit and the cavity, respectively. We denote the product of two qubit- 
cavity states as \qrfn) = 1 77) 1 ® \rf)2 ® \n) C9 where \rf)j label the j-th qubit state (\g) or |e)) and n 
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labels the cavity Fock state. 

The first signature of tripartite interactions is revealed by spectroscopic measurements^!! as a 
function of the detuning Ai^/h = uj\ — uo c of qubit 1 when qubit 2 and the cavity are resonant 
(cj 2 = w c ). In the case of a single qubit-cavity system, the Jaynes-Cummings model predicts a 
single vacuum Rabi-mode splitting of the qubit stated Here, the single qubit states are split twice 
by the mutual interaction of all three systems, as shown in Fig. [T}3. We can interpret this as due 
to the coupling between the bare qubit 1 and the antisymmetric pair of maximally entangled Bell 
states between qubit 2 and the cavity. The two avoided crossings in the spectrum occur along the 
qubit 1 detuning curve, symmetrically displaced about the tripartite resonance (u)\ = uj 2 = cj c ). 
These measured curves agree well with a full analysis of the two-qubit Jaynes-Cummings or Tavis- 
Cumming^3 model. 

With independent control over both qubits, we can easily explore a convenient state-space whereby 
a single photon of energy huj c is shared by our tripartite system. Using a similar technique estab- 
lished for inducing coherent interactions between a single qubit and a cavitjPl we investigate the 
evolution of vacuum Rabi oscillations between qubit 1 and the cavity as a function of the detuning 
of qubit 2 A 2 , c //i = cj 2 — u c from the joint qubit 1 -cavity system {ui = cj c ). For simplicity, we use 
the term "photon" even when describing a single excitation in the qubit. We begin with both qubits 
in their ground state and qubit 1 far off-resonance from the cavity (see pulse diagram in Figj2^), 
then we excite qubit 1 with a photon using a tt pulse and bring it onto resonance with the cavity 
(using a shift pulse) for a given evolution time period t e followed by simultaneous measurement of 
both qubifcP^I. When qubit 2 is far enough off-resonant, the resultant vacuum Rabi oscillations are 
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characterized by the frequency Q = 2 g, as seen on either side of Fig.[2]3,d. Here, the exchange 
between qubit 2 and the qubit 1 -cavity system is energetically prohibited, so that qubit 1 under- 
goes basic vacuum Rabi oscillations with the cavity alone. However, when all three systems are 
on-resonance with each other, the photon begins in qubit 1 and then 'spreads out' to the cavity until 
becoming shared also with qubit 2, eventually moving completely to qubit 2. As time progresses, 
the photon eventually returns to qubit 1. 

In this anti- symmetric mode, the oscillation frequency is given by Q a = Q /y/2. As the system 
evolves, the photon is never completely transferred to the cavity. There are times when the photon 
is entirely in qubit 1 or entirely in qubit 2, otherwise the system occupies a continuum of entangled 
states of both qubits and the cavity (see Fig. (2jT). By measuring the two qubit simultaneously^!, we 
can extract the joint probabilities P eg0 and P ge0 for single-photon states \egO) and \geO), respec- 
tively. A theoretical model including the finite rise-time of the shift pulse (~ 10 ns) agrees well 
with the experimental data (Fig.[2jl-f). 

The above experiment lends itself to a simple geometric description that can help us visualize the 
system dynamics. By use of equation (1), we can identify the unitary evolution U (t) = e~ lHt / h of 
the system with a three-dimensional rotation R n (tp) = e~ m ' X( ^ about n = (0, g 2 , —gi)/ \fg\ + 9% 
with ip = \/g\ + g\ t and X = (Xi, X 2: X 3 ). Here, (X k )ij = — ie^ helps generate the rotation, 
and €ij k is the totally antisymmetric Levi-Civita tensor. Time evolution of the system then corre- 
sponds to orbits on a unit sphere azimuthal to the vector n, where (x, y, z) (\ggl), \eg0), \ge0)) 9 
as shown in Fig.|2^,h. By taking the amplitudes of the three coupled states as real, absorbing any 
overall phase into a redefinition of the states, we can construct a (unit) state vector analogous to 
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that used for a single spin- 1/2 system on the Bloch sphere. In this case, as the state vector precesses 
about n and away from any of the coordinate axes, entanglement evolves over time between all 
three systems. For the experiment described above, we start with an initial condition correspond- 
ing to the state \egO). When qubit 2 is far off-resonance (Fig.|2^), the system precesses at Q Q about 
n = (0, 0, 1), showing simple vacuum Rabi oscillations between qubit 1 and the cavity involving 
the states \eg0) and \ggl), generating bipartite entanglement. However, when all three systems 
are on-resonance (uoi = uj 2 = lj c ), H g 2 X 2 giX 3 g{X 2 — X 3 ), n (0, 1, — 1)/\/2, and 
(p = V2gt leading to a "tripartite evolution". Now the initial state vector \eg0) precesses about 
n so that the trajectory passes from the |e#0)-axis into a region where the photon is shared with 
the cavity and then through the — |ge0)-axis (see Fig.[2]i). The oscillations in the two qubits then 
follow the anti- symmetric mode frequency fi a . This single-photon "tripartite sphere" representa- 
tion provides an intuitive picture for visualizing the equivalence between entangled states in the 
same classP^In this case, the local operations are vacuum Rabi oscillations or tripartite evolution. 
We can see that any arbitrary single-photon tripartite state can be created and subsequently trans- 
formed into any other state on the tripartite sphere, much like unitary operations and rotations on 
the Bloch sphere. Of particular interest is the fact that a specific initial state will follow a specific 
trajectory under tripartite evolution, transforming the amount of entanglement continuously. Be- 
low, we determine the conditions for directly demonstrating transformations between Bell and W 
states, starting from an initially pure state. 

We begin with a single photon in qubit 1 or qubit 2. As shown above, vacuum Rabi oscillations 
represent arbitrary rotations in the \ggl)-\eg0) plane (between the cavity and qubit 1) or the \ggl)- 
\ge0) plane (cavity and qubit 2). These two operations in succession allow us complete access 
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to the \egO)-\geO) plane, and, thus, the ability to prepare any initial state on the entire single- 
photon tripartite sphere. In order to generate Bell and W states, we can start with the photon in 
the cavity, \ggl). Under tripartite evolution the system passes first through the W state, |W) = 
(-|ggl) + |egO) + |ge0))/\/3, and then through the Bell state, |Bell) = -(|egO) + |ge0))/\/2, as 
the system vector rotates about the n- vector, n = (0, 1, — l)/y/2 as shown in Fig.^. In total, the 
system will pass through two Bell states and four W states for one full revolution about n. The 
frequency fl s = \f2 Q of qubit oscillations follows from the definition of and the arc traced 
out by the system trajectory. In this symmetric mode (Q s = 2Q a ) the photon "splits" as it leaves 
the cavity, having an equal probability for going to qubit 1 or qubit 2, and subsequently returning 
completely to the cavity. 

Experimentally, we sample a variety of initial states by allowing qubit 1 (which initially has the 
photon) to undergo vacuum Rabi oscillations with the cavity for a delay time period t d before we 
bring qubit 2 into tripartite resonance. Fig. [3]>e shows a prediction for the unitary evolution of the 
system for nearly a continuum of values for i d . Here, the joint probabilities are P gg i 9 P eg o, and 
P ge0 for states \ggl), \eg0), and \ge0), respectively. Notice that for t d = 2tt/Qo, the system will 
exhibit the anti- symmetric mode (indicated along the dashed line) as described earlier. However 
when t d = tt/Qo 9 we prepare the (initial condition) \ggl), allowing for a tripartite evolution of the 
symmetric mode. After a period of time t e = 7r/4^ s , the excited-state probability for both qubits 
is 1/3 and the system is in the | W) -state, with the photon equally distributed among the two qubits 
and cavity. After a period of time t e = n/2Q s , the excited- state probability for both qubits is 1/2 
and the system is in the bipartite | Bell) -state. These points are indicated in Fig.[3]:-e, with the first 
three states shown as vectors on the tripartite sphere in Fig. [3j}. Here the simulations have included 
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finite energy relaxation and the rise-time of the shift pulses. 

We simultaneously measure both qubits and observe the occupation probabilities of the two qubits 
over time as they evolve from a continuum of initial states, superposition states of qubit 1 and the 
cavity. Although possible, as explained later, we do not measure the cavity state directly. Fig.|4j:,d 
shows extracted line cuts from Fig.|4^,b for two initial conditions (dashed lines) corresponding 
to the symmetric and anti-symmetric modes. As can been seen from Fig.|3]l,e, the theoretical 
predictions for the evolutions agree with the measurements. For the symmetric mode we find in- 
phase oscillations of the two qubits at £} s ~ Q , while for the anti- symmetric mode, we find 
that the two qubits oscillate out of phase with each other with the anti- symmetric mode frequency 
Q a ~ Q Q / y/2, where Q is the frequency of the vacuum Rabi oscillations that occur during the 
delay time period t d (lower right hand corner of Fig. [4^). The measured frequencies agree within 
~ 15 % of the ideal case, due to the finite rise-time of the shift pulses and some residual nonzero 
detuning of each qubit frequency. A theoretical model including these imperfections (solid lines) 
agrees well with the data. 

In the present experiment, we improved the previous design 14 27 by reducing the qubit junction 
areas to reduce the number of two-level system defects. This more than doubled the qubit visibility 
and provided the necessary 'clean' cavity region for observing tripartite interactions. However, it 
was not possible to perform two-qubit state tomography over the required time scales due to short 
relaxation timeP^ matched with the continued presence of two-level system defects that limited 
the qubit visibility to < 50 %. With further reductions in junction size, we can raise the single 
qubit visibility to 90% allowing full tomographic characterization of both qubitsP! 
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In the future, we intend to perform correlated measurements and tomography of this tripartite 
system. This requires a fast, dispersive measurement and readout of the qubits to solve three 
difficulties. First, the tunneling-based measurement of either qubit will populate the cavity with 
unwanted photons due to a crosstalk processpl Second, a dispersive measurement will increase 
qubit visibility ensuring clear tomography. And third, after measurement of the two qubits, subse- 
quent qubit rotations will ensure proper state preparation for one of the qubits, making it ready for 
re-interaction with the cavity. In this way, we can reuse one of the qubits through state transfers', 
to fully determine the cavity state 30 . Improvements are currently underway to modify our slow 
switching-current SQUID readout to a fast, dispersive resonant readout. 

Tripartite interactions provide a means of engineering qubit entangled states. We have provided a 
theoretical description of the two-qubit Jaynes-Cummings or Tavis-Cummings model and a free- 
evolution protocol for the deterministic preparation of Bell, W, or arbitrary entangled states with 
one shared photon. This includes a convenient geometric description helpful for visualizing the 
unitary evolution. We have presented clear experimental evidence of tripartite interactions between 
two 
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superconducting phase qubits and a resonant microwave cavity. Our system has a fundamental 
advantage over strictly multi-qubit systemP^^l The cavity provides a larger degree of freedom as 
opposed to the two discrete levels associated with qubits and can be thought of as a multiphoton 
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register or entanglement resource. Arbitrary preparation of multiphoton states in this cavity via 
one of the qubits 30 , and subsequent interactions for entanglement distribution will allow for the 
creation of further classes of entanglement, such as a GHZ state. 
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Figure 1 : Circuit and spectroscopy, a, Optical micrograph of the electrical circuit with two Joseph- 
son phase qubits (qubit 1 inset overlay right), each with loop inductance ^700 pH and critical cur- 
rent ~ 0.91 /iA (junction areas ~ 6 /xm 2 ) shunted by use of interdigitated capacitors (Cj ~ 0.7 pF, 
including junction capacitance) with vacuum gap crossovers (inset overlay left), capacitively cou- 
pled (C c ~ 6.2 fF) to a coplanar waveguide resonant cavity (of full length ~ 7 mm). The device 
was fabricated with standard optical lithography with A1/A10 X /A1 junctions on a sapphire sub- 
strate, with Si0 2 as an insulator surrounding the junctions, b, Microwave spectroscopy of qubit 1 
as a function of detuning Ai jC = uj\ — uj c with uj 2 = u) c - Ai jC is varied through the d.c. flux bias 
coils and qubit 1 is excited by microwaves applied through the m.w. (microwave) coil (seen in a). 
The intensity color scale represents the probability of qubit 1 tunneling after the measure pulse. 
The dashed diagonal line shows the bare qubit 1 transition frequency. The dashed horizontal lines 
represent the two maximally entangled Bell states between qubit 2 and the cavity. 
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Figure 2: Demonstration of basic tripartite interactions, a, Description for creating a photon in 
qubit 1 by use of a tt pulse, then shifting (solid line) onto resonance with the cavity and qubit 2 for 
various qubit 2 detunings (dashed line). After an evolution time period t e the qubits are measured 
simultaneously. b,c, Measured excited state joint probabilities P eg0 and P geQ for states \egO) and 
\geO) 9 respectively, during tripartite interactions after qubit 1 has been excited by a tt pulse and 
shifted onto resonance with the cavity as a function of the detuning A 2 , c = uj 2 — u c of qubit 2. 
d,e, Theoretical predictions including qubit visibility and finite rise-time of the shift pulse. The 
resulting asymmetry is relatively insensitive to the exact shape of the shift pulse and can be at- 
tributed to additional interference due to finite detuning of ~ 0.3#. Here, we used an exponential 
rise-time of ~ 10 ns. f, Line cut of the on resonance tripartite interactions with corresponding the- 
oretical prediction (solid line). g,h, The red arrow provides a visual cue to the circular trajectory 
of the tripartite vector, g, Tripartite sphere representation during simple vacuum Rabi oscillations 
of qubit 1. h, Tripartite sphere representation during the tripartite evolution from the initial state 
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Figure 3: Experimental protocol and theoretical predictions for generating generalized arbitrary 
single photon tripartite evolutions, a, Pulse sequence: (1) A photon is inserted in the system 
by exciting qubit 1. (2) A shift pulse brings qubit 1 onto resonance with the cavity, producing 
vacuum Rabi oscillations. (3) A shift pulse brings qubit 2 onto resonance after the delay time £ d , 
initiating tripartite interactions that evolve over a time period t e — t d . (4) Both qubits are measured 
simultaneously, b, Tripartite sphere representation of the tripartite evolution for the initial state 
\ggl) prepared during a delay time period t d = n/Q . The red arrow provides a visual cue to 
the circular trajectory of the tripartite vector, c, Predicted state occupation of one photon in the 
resonant cavity. d,e, Predicted joint state probabilities P eg0 and P ge0 for measurement of qubit 1 
and 2 as functions of both t d and t e . Blue color represents low values, red represents high. The 
simulations were performed with a finite energy relaxation time, a finite shift pulse rise-time (~ 10 
ns), and some residual qubit detuning of ~ 0.3g from the cavity. 
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Figure 4: Experimental demonstration of arbitrary tripartite interactions between both phase qubits 
and the cavity. a,b, Measured joint state probabilities P eg0 and P ge0 for measurement of qubit 1 and 
2 as functions of both t d and t e . c, Extracted curves (along dashed line) for the initial state \ggl) 
producing a tripartite evolution of the symmetric mode, showing in phase oscillations along with 
theoretical predictions (solid lines) from Fig.[3]i,e. During this evolution, the system's entangle- 
ment continuously transforms starting from a pure state \ggl). In the ideal case, the system evolves 
through W and Bell states, d, Extracted curves (along dashed line) for the initial state \egO) pro- 
ducing a tripartite evolution of the anti- symmetric mode, showing out of phase oscillations along 
with theoretical predictions (solid lines) from Fig.|3jl,e. 
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